Calculation of Some Expected Values for Parameterized Mean 

Model with Gaussian Noise 
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Consider the measurement model 
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y = g( x ) + v 



(i) 



where 



• i £ Mis the unknown scalar we would like to estimate described by the prior distribution Af(x; 0; a 2 ) 
where the notation • ; x, E x ) denotes a (real-variate) Gaussian density with mean x and covari- 
ance E a . 

• y £ M™ y is, in general, a complex measurement vector; 

• g( ■ ) : R — > C"" is, in general, a complex-valued observation function; 

• d£ C" m is circular symmetric complex Gaussian measurement noise with zero-mean and covariance 
<JyIn y where I n denotes an identity matrix of size n y x n v . The noise v is assumed independent of 



In this document we are going to derive analytical formulae for the following functions 
Lh,,x(si,s 2 ,hi,h 2 ) —E mm(L Sl (y,x + h u x),l)mh\(L S2 (y,x + h 2 ,x),l) , 
Vy\x(si,s 2 , hi,h 2 ,x) =E y min(Lf (y, x + hi,x),l) mm(£f (y, x + h 2 ,x),l) 
Hx{si,s 2 ,hi,h 2 ) =E X mm(L 2 !1 (x + hi,x),l) iam{L s 2 :! {x + h 2 ,x),l) 
where si,s 2 £ Z, hi, hi € R™* and 



,Py, x (y,x) A Py\x{y\x) A p x {x) 



(2) 
(3) 
(4) 

(5) 



1 Calculation of /i y;a ;(si, S2, /ii, /i2) 

We first calculate L(y, x + h,x) as 

CM{y; g(x + h),o 2 v I ny )M{x + h; 0; a 2 x ) 



L(y, x + h, x) 



CAf{y;g(x),a 2 I ny )Af(x;0;a 2 
2 

:exp — 1Z \v R d(x, h)\ — b(x, h) 



(6) 
(7) 



where the notation CAf(y; y, H y ) denotes a circular symmetric complex Gaussian density with mean y 



and covariance E y and 



d(x,h) =g(x + h)- g(x), 

b(x,h) ^\\d(x,h)\\ 2 + \x T h+ ^\\h\\\ 



(8) 
(9) 
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1 



Then we have 



/%,z(si,s 2 , hi, ha) 
min exp 



— ^-1Z {v H d(x, hi)} — sib(x, hi) 

, 1 )p(v) dvp(x) dx. 



x mm exp 



—^-TZ {v H d(x, h 2 )} — s 2 b{x, h 2 ) 



(10) 



We are first going to handle the inner integral on the right hand side of (|10l) which we call as Zi as 
follows. 



X\ = J min I exp 
x min I exp 



^K{v H d(x,hi)} - Sl b(x,hi) ,1] 

1 )p(v) dv 



-^-TZ{v H d(x, h 2 )} — s 2 b(x, h 2 ) 



p(v) dv 



Vi 



exp 



v 2 



2s 2 



K{v K d(x, h 2 )} — s 2 b(x, h 2 ) 



exp 



v 3 



exp 



Vi 



— ^TZ{v R d{x, hi)} — sib(x, hi) 
— ^-72.{w H (i(x, hi)} — sib(x, hi) 



p(v) dv 

p(v) dv 
2s 2 

exp 



TZ{v a d{x, h 2 )} - s 2 b(x, h 2 ) 



where the sets Vi, V 2 , V3 and V4 are defined as follows. 
Vi = £ C 
V 2 ^[v£ C 
V 3 = {v G C 

V 4 = {ve C n " n{v^d(xM)}< alH l M ^n{v^d(xM)}< ' iH l' h2) y 
Substituting (0) and the identity p(v) — CN{v; 0, o 2 I ny ) into (TTT1) . we get 
Zi =P(«G Vi\v~CN(v;0,al)) 





h 2 ) 


2 






h 2 ) 


2 






h 2 ) 


2 






h 2 ) 



p(v) dv (11) 

(12a) 
(12b) 
(12c) 
(12d) 



+ exp 
+ exp 
+ exp 



h 2 )\\ 2 - s 2 b{x, h 2 ) 



-4||d(x, hi)\\ 2 - sib(x, hi) 



P(veV 2 \v~ CN{v- s 2 d(x, h 2 ), al)) 
P(v e V 3 \v ~ CN(v; sid(x, hi), a 2 v )) 



■\\sid(x, hi) + s 2 d(x, h 2 )\\ 2 - sib(x, hi) - s 2 b(x, h 2 ) 



x P(y e Vi\v ~ CAf(v; sid(x, hi) + s 2 d(x,h 2 ),al)). 



(13) 



We now define the real scalars ai and a 2 as ai = 1Z {v li d(x 7 hi)} and a 2 = 1Z {u H d(x, h 2 )} . Since each 
of the probabilities on the right hand side of (|13|) are conditioned on v being distributed with a circular 
symmetric complex Gaussian density and since ai, a 2 are linearly dependent on v, we have the vector 
a = [ai,a 2 ] T distributed with a Gaussian density which gives 



h =P 



ai> - 
v<*2>- 



2 

rgh(g,h 2 ) 



a~ N [ a;a[(x,hi,h 2 ), -^-T(x,hi,h 2 ) 



2 



+ exp 
+ cxp 
+ exp 
x P 



h 2 )\\ 2 - s 2 b(x, h 2 ) 



1 

Ql < 
a 2 < 



\\\d(xM)\\ 2 - sib(xM) 



ai > 



\a 2 < — 



ai< 



a 2 ^ o 



(/ ~ JV ( a; 03(2;, hi, h 2 ), -rr T ( x , hi,h 2 ) 



a ~ M ( a; 03(0:, /ii, h 2 ), "y r ( x ' M 



||sid(x, /ii) + s 2 d(x, h 2 )\\ 2 — sib(x, hi) — s 2 b(x, h 2 ) 



2 

rgt-(x,h 2 ) 



a ~ AM a; a' 4 (x, hi,h 2 ), -£-T(x, hi,h 2 ) 



(14) 



where 

a' 1 (x,hi,h2) = [0,0] T , 



a' 2 (x, hi, h 2 ) =s 2 \jZ{d n {x, hi)d(x, h 2 )}, \\d(x, h 2 )\\ 



(15a) 
(15b) 

a' 3 (x,hi,h 2 ) =si [\\d(x,hi)f,K{d Ji (x,hi)d(x,h2)}] X , (15c) 

5i(x, hi, h 2 ) 4 [ Sl ||d(x, /ii)|| 2 + s 2 ft{d H (a;, hi)d(x, /i 2 )}, *a||d(ss, Mil 2 + si^{d H (x, ftijdfo ^ 2 )}] T , 

(15d) 



T(x,h 1 ,h 2 ) = 



||d(a:,MI| 2 n{d^{x,hi)d{x,h 2 )} 
K{d H (x,hi)d(x,h 2 )} \\d(x,h 2 )\\ 2 



(16) 



Each of the probabilities on the right hand side of (|14[) can be written using the cumulative distribution 
function A/cdf 2 ( •,-,•) of a bivariate Gaussian random variable to give 



h =/Vcdf 2 
+ cxp 

x yvcdf 2 



a 2 b(x,hi) a 2 b(x,h 2 ) 



:\\d(x, h 2 )\\ 2 - s 2 b(x, h 2 ) 
a 2 v b(x,hi) a 2 b(x,h 2 ) 



+ exp 
x 7Vcdf 2 



-±\\d(x,hi)\\ 2 - sib{xM) 



r cr 2 6(a;, hi) a 2 v b(x,h 2 ) 



; ai(x, hi, h 2 ), -y-r(x, hi,h 2 ) 



a 2 

; a 2 (x, hi,h 2 ), yr(a:, hi,h 2 ) 



; 03(2;, hi, hi), -^ r (a;, hi, h 2 ) 



+ exp 
x 7Vcdf 2 



—?\\sid(x, hi) + s 2 d(x, h 2 )\\ 2 - sibtx, hi) - s 2 b(x, h 2 ) 
at 



a 2 b(x,h x ) a 2 b(x,h 2 ) 



; a 4 (x, hi,h 2 ), yr(i, fc^, /i 2 ) 



where 



ai(x, /ii,/i 2 ) =-ai(a;,/ix,/i 2 ), 
a 2 (x,hi,h 2 ) = ^ ° d' 2 (x,hi,h 2 ), 

da(x, hi,h 2 ) = _ x d' 3 (x,hi,h 2 ), 
64 (x, hi,h 2 ) =d' 4 (x, hi,h 2 ), 



T(x, hi, ha) = 



\\d(x,hi)\\ 2 -n{d n (x,hi)d(x,h 2 )} 
-K{d H (x, h x )d(x, h 2 )} \\d(x, h 2 )\\ 2 



(17) 

(18a) 
(18b) 

(18c) 
(18d) 

(19) 



3 



Using © in ([T7|). we get 
h =A/"cdf 2 



crgb(a:, fei) ag&(cc, h 2 ) n T 
2 ' 2 



ai(a;, hi,h 2 ), ~f r ( x i h u h 2 ) 



+ exp 
x A/cdf 2 



S 2 _ s 2 n ,/ , x M 2 S 2 Tl S 2 m, ,|2 



cr^6(x,/ii) a 2 b{x,h 2 ) 1 T _ , , , x ^ 
;a 2 (a:,/ii,ft 2 ), —T(x,hi,h 2 ) 



+ exp 
x A/cdf 2 
+ exp 



2 ' 2 
alb(x,hi) <Jyb(x, h 2 ) 



n r 2<j2 

T ^ 2 

; a 3 (a;, hi,h 2 ), -yT(x, /ii, /i 2 ) 



2 

2sis 2 . 



Sl/tl + S 2 ft 2 Si ,,, ||2 S 2 2 



x A/cdf 2 
=A/cdf 2 



2 ' 2 



a%b(x,hi) a 2 u b(x,h 2 ) 



+ exp 

x A/cdf 2 

+ exp 

x A/cdf 2 

+ exp 
x exp 

x A/cdf 2 



s| - s 2 



\\d{xM)\? + S -^\\h 2 f 



; a 4 (a;, h 1} h 2 ) 7 -ylXx, /ii, ft 2 ) 



0i(x, hi,h 2 ), -fT(x, hi,h 2 ) 
exp 



.2- ..a-^jIlM 



<r 2 b(x, ft x ) ggb(a;, fe 2 ) 
2 ' 2 



; a 2 (x, hi,h 2 ), hi,h 2 ) 



2t 2 

a%b{x,hi) olb(x,h 2 y T 



s 2 - Si 



2 

„2 



exp 

a s (x, hi,h 2 ), -^-^(x, hi,h 2 ) 
2sis 2 , 



., \\d(x, h^f + ?*-Jl\\d(x, h 2 )\\ 2 + -^TZ {d H (x, h 1 )d{x, h 2 )} 

(7 G 

V V V 



£^ l|Alf + £^ IM > + -, 



■hjh 2 



a%b(x,hi) cr 2 v b(x, h 2 ) 
2 ' 2 

Substituting the result (|2~TT) into (ITU1) and 



exp 
a 2 



(20) 



si/ii + s 2 h 2 1 

5 x - —j si/ii + s 2 h 2 \ 



; a 4 (x, hi,h 2 ), -^-^(x, hi,h 2 ) 



(21) 









^y,x(si,S 2 ,hl,h 2 ) =£'a^(x;0,o-2) 


A/cdf 2 ^ 


2 

cr^ b{x.h2) 
L 2 J 



carrying out straightforward algebra, we obtain 
ai(x, hi,h 2 ), r(x, ftx, h 2 ) 



x A/cdf 2 



exp 



S 2 ~ S 2 || ,/ , x |, 2 , S 2 - s 2 

— \\d(x,h 2 )\\ 



^-i,M 2 



~~ 2 

<7 v b(xjl2) 

2 



; a 2 (x, hx,h 2 ), ~f^( x , hi,h 2 ) 



4 



+ E Af(x;- Sl h 1 ,a^) 

x TVcdf 2 



exp 



at 



-\\d{xM)\\ 2 + 



2al 



-\\hi 



a v b(x,h2) 
2 



^2 

; a 3 (x, h 1 ,h 2 ), -y r ( x ' 



+ £ jV(s;-(si/n+s 2 /i2),<^) 

2sis 2 



exp 



+ 



-^{d H (x,/ii)d(x,/i 2 )} 



exp 



S l — S l II z, \\2 , S 2 _ s 2 || , || 2 SiS 2 



2ct? 



-lift 



2a? 



ft 2 r + 



x TVcdf 2 



cr^fc(a:,/ii) 
2 

2 



; di(x, hi,h 2 ), Y r ^' 



2 Calculation of fi y \ x (si, S2, hi, x) and /i x (si, S2, /ii, ^2) 

We first calculate Li( • , • , • ) and L 2 ( ■ , • ) as follows. 

-j,\\y-g(x + h)t 



Li{y,x + h,x) 



exp 



L 2 (x + h, x) 



exp 
= exp 

= exp 

= exp 
exp 



-TZ {y^d{ Xl h)} - -\\g(x + h)f + -\\g{x) 



K{v R d{x,h)} - —\\d(x,h)\\ 



n{v R d{x,h)} - h(x,h) 



exp 



= exp 



=x L h- 



2aV 



exp [-b 2 (x,h)] , 



where 



h(x,h) =^\\d(x,h)\\ 2 , 



Then the function fj, y \ x ( ■ , • , • , • , • ) is given as 
fJ- v \x{si,s 2 ,hi,h 2 ,x) 

= J min ^cxp —^-H {v H d(x, hi)] — sibi(x, hi) , l^j 



x mm [ exp 
dv 



^-IZ {v H d(x, h 2 )} - s 2 bi(x, h 2 ) 



, 1 ) p(v) dv 



= / P(v) 

JVi 



exp 



v 2 



2s 2 



TZ {v H d(x, h 2 )} - s 2 bi(x, h 2 ) 



p(v) dv 



/ exp 

Jv 3 



-^-7?. {v H d(x, hi)} — s\b\{x, hi) 



exp 



Vi 



2s x , 



p(v) dv 
~2s 2 

exp 



-71 {v H d(x, hi)} — sibi(x, hi) 

where the sets Vi, V 2 , V 3 and V4 are defined as follows. 

Vi = (v e C n » n{v"d(x,h 1 )}> ^ bl( 2 thl) &K{v ll d(x,h 2 )} 



1Z {v H d(x, h 2 )} - s 2 bi(x, h 2 ) 



V 2 ^ 



}■ 

ju e o K{t, H rf( a; ,h 1 )}> g "'' i S ; 3! '' ii) fcTC{ t) H rf( a; ,/i 2 )}< g "'' i( 2 3; -'' 2) |, 





/, 2 > 


2 






^2) 


2 






h 2 ) 


2 






h 2 ) 



p(v) dv (32) 

(33a) 
(33b) 
(33c) 
(33d) 



Substituting ((29]) and the identity p(v) = CN(v; 0, a 2 I ny ) into ([32]). we get 



fJ- v \x(si,s 2 ,hi,h 2 ,x) 
=P 



K{v H d(x,hi)}> °* H 2 * M 



n{v n d(x,h 2 )}: 



?ji>i(x,h 2 ) 



v ~ CN{v;Q,all nv ) 



+ exp 

+ exp 
+ exp 
x P 



V 4 ~ S 2 



r s\ - Si 



\\d(x,h 2 )\\ 



\\d(x,hi)\\ 



p 



p\ 



si - Si 



TZ{v H d(xM)}l 



TC{jj H d(x,ft, 2 )}; 

s| - s 2 







2 






<»2> 


2 








2 






h 2 ) 



v ~ CAf{v; s 2 d(x,h 2 ),a 2 I ny ) 
v ~ CAr(v;sid(x,hi),<jll ny ) 



\\d(x, hi)\\ 2 + ?»-Jl\\ d ( x , h 2 )f + ^lp-U{d H (x, hi)d{x, h 2 )} 



at 



Tl{v a d{x,h 1 )\' 



K K{v"d(x,h 2 )}< a " bl{ 2 th2) 

Continuing in the same way as in Section [TJ we get 



v ~ CAf(v; sid(x, hi) + s 2 d(x, h 2 ),all ny ) 



^ y \x(si,s 2 , hi, h 2 , x) =A/cdf 2 



cr y bi(x,hi) 
~~ 2 
c y b 1 (x,h 2 ) 
2 



; ai{x, hi,h 2 ), -£T(x, hi, h 2 ) 



+ exp 
x A/cdf 2 



4 ~ S 2 



+ exp 
x A/cdf 2 



s{ - Sl 



\\d(x,h 2 )\\ 

cr^fci (x,hi) 

~~ 2 

<J v bi(x.h 2 ) 

2 

\\d(x, hi)\\ 



& "tbl (x,hi) 
2 

cr^&i (x,h 2 ) 



; a 2 (x, hi,h 2 ), ~^^{x, h x ,h 2 ) 



; a 3 (x, hi,h 2 ), -y r (a;, hi, h 2 ) 



+ exp 
x 7Vcdf 2 



T s\ - Si 



(34) 



(35) 



\\d(x, hi)\ 



! 2 + $-^-\\d(x, h 2 )\\ 2 + {d H (x, hi)d(x, h 2 )} 



<J y bi(x.hi) 
2 

<y v bx{x,h 2 ) 
2 



; a 4 (x, hi,h 2 ), -^-"^(x, hi, h 2 ) 



(36) 



G 



Similarly the function fx x (•,-,•,•) is given as 

Mx(si,S 2 ,/ll,/l2) 

min (exp [— si^x, ^-l)] , 1) rnin (exp [— s 2 b 2 (x, h 2 )] , l)p(x) dx 



(37) 



p(x) dx + , exp \— sobiix, h 2 )] p(x) dx 

b 2 (xM)<0 i y ' Jba(x,hi)<0 1 K n K ' 

b 2 (x,h 2 )<0 b 2 (x,h 2 )>0 

exp [— sib 2 (x, h{)]p(x) dx + 



b 2 (x,h 1 )>0 
b 2 (x,h 2 )<0 



b 2 (x,hi)>0 
b 2 (x,h 2 )>0 



exp {-sib 2 (x, hi)] exp [-s 2 b 2 (x, h 2 )]p{x) dx 



(38) 



--P 



(b 2 {x,h 1 )<0 
\b 2 (x,h 2 )<0 



+ exp 
+ exp 
+ exp 
--P 



2ol 



2ol 



x ~ N(x]0;al)j 

/ b 2 (x,h! 
\b 2 (x,h 2 



n 



i)>0 

•,ft 2 )<o 



S l ~ Sl II A 112 , S 2 - s 2 || , „2 . s l s 2 , T , 
"ill ^T^ - IIM H ^~h 1 h 2 



2ol 

CLl < 



A/"(x; — si/ii, (T x ) 



)>0 
h 2 )>0 



x ~ A/"(x; — si/ii — S2/12, o"x) ) (39) 



a 2 < 
+ exp 



TV a;a'i(hi,h 2 ), -^A(h 1 ,h 2 ) 



+ exp 
+ exp 
x P\ 



2a? 



Sl ~ Sl »hi\\ 2 



2ol 



PI 



P 



ai < 
a 2 > 

ai > 
a 2 < 



a — A/" ( a; a' 2 (hi,h 2 ), —~A(hi,h 2 ) 



a ~ ( a; a^/ii, /i 2 ), — ^ A(/ii, /12) 



S l _Sl IIA 112 , S| — S2 |. 2 «1S2,T. 

"ill H ^-5— 11^2 1 1 H 2~h 1 h 2 



ai > 
a 2 > 



2a? 



a ~ AM a; a' 4 (hi,h 2 ), —A(hi,h 2 ) 



where 



a' 2 (/n,M [-2s 2 ^/i2 + ||/M|| 2 , (1 - 2,s 2 )||/i2|| 2 ] r 



a' 3 (h u h 2 ) [(l-2 Sl )||M 2 ,-2 Sl ^ 2 + ||M2] T : 



a 4 (hi,h 2 ) 



1 

; 2^! 



(l-2si)||/li|| 2 -2s 2 /l^2 

-2sihjh 2 + {l-2s 2 )\\h 2 \\ 2 



(40) 

(41a) 
(41b) 
(41c) 
(41d) 



A(hi,h 2 



\\hi\\ 2 hjh 2 
hjh 2 \\h 2 f 



(42) 



Each of the probabilities on the right hand side of (|40j) can be written using the cumulative distribution 
function Afcd{ 2 ( ■,-,-) as 



Vx{si,s 2 ,hi,h 2 ) 



=A/bdf 2 





' " 








;ai(hi,h 2 ), —K{hi,h 2 ) 



7 



+ exp 
+ exp 
+ exp 



4 ~ S2 | 

2a? 1 



s l - s l 

2a? 



h 2 \\ 2 



hi\ 



A/"cdf 2 
7Vcdf 2 



;a 2 (fti,ft 2 ),— A(ft 1; ft 2 ) 



■,a 3 (hi,h 2 ),\A(hi,h 2 ) 



*N* l|IM 2 + £|zi2| N | a + ^T^ 



2a? 



2a? 



A/bdfo 



where 



ai(fti,ft 2 ) =ai(/ii, ft 2 ), 

2 2 (/ll,/l2)= o -1 5 2 (ftl,^2), 

a 3 (hi,h 2 )= ~ x ° a' 3 (ln,h 2 ), 
a,4{hi,h 2 ) = - a' 4 (hi,h2), 



A(h u h 2 



|fti|| 2 -hjh 2 
-hjh 2 \\h 2 \\ 2 



i 04(^1,^2), -J A(fti, /is) (43) 
ct 2 / 

(44a) 
(44b) 

(44c) 
(44d) 

(45) 



3 Special Cases 



In this section, we are going to find the expressions for the following quantities one by one: fi y ^ x (l, 1, ft, ft), 
fi y , x (l,0,h,h), fi y \ x (l,l,h,h,x), fi y \ x (l,0,h,h,x), fi x (l,l,h,h), ^(1,0, ft, ft). 

• M3/,x(li 1: M : Substituting sx = s 2 = 1 and fti = h 2 = h in (|22[) . we get 

CT 2 N 
;a x ' (a, ft, ft), -£T(x,h,h) 







a1b(x.h) 


Hy,x(l, 1, h, ft) =-E\/'(x;0,<7g) 


■A/cdf 2 ^ 


2 

2 J 











A/cdfa ^ 


2 

cr 2 b(x,h) 
L 2 J 






<j 2 v b(x,h) 


E M{x;-h,al) 


A/bdf 2 ^ 


2 

L 2 J 



;a2' x (a;,/i,/!,), -^-Y(x,h,h) 



;al A (x,h,h),Y T ( x , h , h ) 



+ ^JV(k;-2/i,ctJ 

x A/cdf 2 



exp 



exp 



a 2 b(x,h) 
2 

a 2 v b(x,h) 
2 



1 (as,/i, ft), yr(a),/i,/i) 



where 



aj :1 (ir, ft, ft) = [0, 0] T , 
a2 4 (a;,/i,/i) = [-||d(a;,/i)|| 2 , ||d(x, ft)|| 2 
a 3 ,:l (x, ft, ft) = [||d(x, ft)|| 2 , h)\ 
a l /{x,h,h) =2 [||d(a;,ft)|| 2 ,||d(a;,ft)| 



T 

21 T 



and 



F(x, ft, ft) = 
r(x, ft, ft) = 



||d(a;,ft)|| 2 \\d(x,h)\\ 2 
\\d(x,h)\\ 2 ||d(a;,ft)|| 2 

||d(x,ft)|| 2 -\\d(x,h)\ 
-\\d{x,h)\\ 2 ||d(x,ft)|| 2 



(46) 

(47a) 
(47b) 
(47c) 
(47d) 

(48a) 
(48b) 



Noting that T(x, h, h) and T(x, h, h) are singular, the quantities related to the bivariate cumulative 
distribution function A/cdf 2 ( • ) can be written as 



A/"cdf 2 y 

A/cdf 2 ^ 
A/bdf 2 ^ 
Afcdf 2 



2 

rgb(x,fe) 



2 

2 

alb(x,h) 
2 

of,b(x,h) 
2 

2 

a 2 v b(x,h) 



; 5i (a;, /i, ft), -^r(x, ft, ft) =Nc&h - 



;a,2{x,h,h),-2-T{x,h,h) \ =0, 



\a3{x,h,h),^-T(x,h,h) \ =0, 



a 2 v b(x,h) a 2 v 



0,^-\\d(x,h)\\ 2 ), (49a) 



(49b) 
(49c) 



; 04(0;, h, h), -^-T(x, h, h) J = A/cdf 1 



(49d) 



where A/cdf 1 ( • ) is the cumulative distribution function for the univariate Gaussian density, which 
gives 



Mj/,a;(l) lj ft, ft) —EfJ( x -fl^l) 



A/cdfi 



ofo(s,ft) n <^ 

2 — ;0,yKx,ft)|| 



exp 



r Hz, /OH 



exp 



rllftll 



A/cdfi 



(T^6(a;, h) 



2\\d{xM\ 2 ^\\d{xM\ 2 



(50) 



fj, y ^ x (l,0,h,h): Substituting si = 1, s 2 = and h\ = li2 = h in l]22p. we get 



/%,x(l, 0, /l, /l) =-EV(x;0,o-J) 



A/cdf 2 



2 

alb(x,h) 
2 







alb(x,h) 


^(a;0,o-J) 


A/cdf 2 ^ 


2 

crjb(x,/l) 
L 2 J 



+ E M(x--hMl) 



Afcdf 2 
Afcdf 2 



2 

al.b(x.h) 
2 

<j 2 v b(x,h) 
2 

a%b(x,h) 
2 



;a{ (x,h,h), -£-T(x,h,h) 



;a 2 ' (x,h,h), -^-T(x,h,h) 



;al' (x,h,h),^-T(x,h,h) 



\aY(x,h,h), -^-Y(x,h,h) 



where 



a\'°(x, ft, h) = [0, 0] T , 
a,2°(x, h, h) = [0, 0] T , 



al'°(x,h,h)±[\\d(x,h)f,-\\d(x,h)\\ 2 ] T , 
dl<°(x,h,h) ^[\\d(x,h)\\ 2 ,\\d(x,h)\\ 2 ] T . 
Again due to the singularity of T(x, ft, h) and T(x, h, h), we have 



(51) 

(52a) 
(52b) 
(52c) 
(52d) 



Mj/,x(l, 0, ft, ft) =£yv( x; o,o2) 



A/cdfi 



-^);0,|||^,MI| 2 



'Af(x;-a,<rg) 



A/cdfi 



cr 2 b(x, h) 



\\d(x,h)\\ 2 ,^\\d(x,h)\\ 



(53) 
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(i y \ x (l, 1, ft, ft, x): Substituting s± = S2 = 1 and h\ = h% = h ii 



we get 



£tj,| x (l, 1, ft, ft, 2;) =A/cdf 2 ( 
+ A/"cdf 2 
+ A/bdf 2 
+ exp 



2 

CT^fel (oi,/l.) 
2 

cr^ 61 (x,h) 
~~ 2 
(7 V bi (x : h) 
. 2 

bi (x,h) 
~2 
a v bi (x,h) 



;a x ' (x,h,h), -±T{x,h,h) 



;a 2 ' (x,h,h), y^/i./i) 



;4' 1 ( a; )^ /l )> y r (a;,/i,ft) 



AAcdf 2 



(T^ 61 (x,h) 
2 

cr^ bi {x,h) 
2 



; 04(0;, ft, ft), -^-r(x, ft, ft) 



Using the singularity of T(x, ft, ft) and T(x, ft, ft), we obtain 

2 ;0, y|M(a:,/i)|| 



^^(l, 1, ft, ft, a;) =A/cdf 1 ^— 



+ exp 



:\\d(x,h)\\ 



AAcdfi 



er 2 &i(a;, ft) 



2|K*,ft)|| 2 ,^|Kz,ft)|| 2 



/Lt y | x (l, 0, ft, ft, x): Substituting Si = 1, s 2 = and hi — ft 2 = ft in (|36p . we get 



/u y | x (l,0,ft,ft,a;) = /Vcdf 2 ^ 
+ 7Vcdf 2 
+ Afcdf 2 
+ A/cdf 2 



(T^bi (x,h) 

<y^b\ (x,h) 
2 

&vbi{x,h) 
~~ 2 
cr v bi(x : h) 
2 

<7^bi(x,h) 
~~2 
<J v b\ (x.h) 
. 2 
cr^ bi (x,h) 
2 

<T V b\ {x.h) 
. 2 



;a\'°(x,h, ft), yr(j,/i, ft) 



;a 2 '°(x, ft,ft), Y r (x,h,h) 



;al'°(x,h,h), -2-T(x,h,h) 



;a4°(x, ft, ft), yr(i,ft,/i) 



Using the singularity of r(x, ft, ft) and T(x, ft, ft), we obtain 



fj, y \ x (l,0,h,h,x) =7Vcdfi 



alh(x,h) a 2 v 



Q,~-f\\d{xM\ 2 



A/bdf 1 



albi{x,h) 



\\d{x,h)\\\^\\d{x,h)\\- 



If we now substitute 61 ( • , • ) from (f29| into (f57| , we get 



^ y | x (l,0, ft, ft, x) =A/cdfi ^~ 
+ A/bdfi 



^^;0,|lM(x,ft)f 



||d(x,ft)|| ; 



Kz,ft)||^|Kx,ft)|| 2 



H*,ft)l 



=2JVcdf a ( - 

' \\ djx, ft) || 

2<7„ 



,0,^||d(z,ft)|| 2 



=1 - erf 



where we used the identity 



ATcdf! (£;£,af) =- 1 + crf 
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jji x (\, 1, ft, ft): Substituting si = S2 = 1 and hi = h% = ft in (j4"3")l . we get 



^(1,1, ft, ft) =7Vcdf 2 

+ A/"cdf 2 
+ JVcdf 2 





' " 








ifi}' 1 ^, ft), -^A(h 1 ,h 2 ) 





' " 


( 







' " 


( 






;al' 1 (h,h),-^A{h,h) 



+ exp 



- 2 M 



A/"cdf 2 



; 04 1 (h, ft), —A(h, ft) 



where 



1 



~a{'\Kh)^— [||ft|| 2 ,||ft|| 2 ] 



-1,1 
a 2 

~i,i 



al>\h,h)±-^[\\h\\ 2 ,-\\h\\ 2 ] T . 



~al'\h,h) ±— [||ft|| 2 ,||ft|| 2 ] 



and 



A(ft,ft) = 
A(ft,ft) = 



\\h\\ 2 \\h\\ 2 
\\h\\ 2 \\hf 

\\h\\ 2 -\\h\ 



-IN 2 IN 2 

Using the singularity of A(ft, ft) and A(x, ft, ft), we obtain 



^(1,1, ft, ft) =A/cdfi 0: 



l|ft|| 2 \\H 



exp 



1 



\\h\\ 



AAcdfi 0; 



3||ft|| 2 \\h\\ 



X " X 



0, ft, ft): Substituting s\ = 1, s 2 = and h\ = /i 2 = ft hi (03]), we get 



^(1,0, ft, ft) = /Vcdf 2 

+ A/bdf 2 

where 





' " 








;a{'°(h,h),^A(h,h) ) + A/"cdi 2 



(62) 

(63a) 
(63b) 
(63c) 
(63d) 

(64a) 
(64b) 

(65) 



;a 2 '°(ft>ft)> ^2 A (ft)ft) 





' " 








= 1,0- - 1-r 



,a3' u (ft,ft),— A(ft,ft) ) + A/"cdi 2 



;al'°(h,h), -^A(h,h) 
at 



-1,0 

a 



al'°(h,h) 
al'°(h,h) 



1 



-||ft|lM|ft|| 2 ] T : 



2^ 2 

a 4 '°(ft,ft) [||ft||M|ft| 



[l|ft|| 2 ,-||ft|| 2 ] T ' 



(66) 

(67a) 
(67b) 
(67c) 
(67d) 
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Using the singularity of A(h, h) and A(x, h, h), we obtain 



=i-«tfjp_ 



(70) 
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